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In this paper, we generalize the ordinary two-photon Jaynes-Cummings model (TPJCM) by con- 
sidering the atom (or ion) to be trapped in a simple harmonic well. A typical setup would be an 
optical cavity containing a single ion in a Paul trap. Due to the inclusion of atomic vibrational 
motion, the atom-field coupling becomes highly nonlinear what brings out quite different behaviors 
for the system dynamics when compared to the ordinary TPJCM. In particular, we derive an effec- 
tive two-photon Hamiltonian with dependence on the number operator of the ion's center-of-mass 
motion. This dependence occurs both in the cavity induced Stark-shifs and in the ion-field cou- 
pling, and its role in the dynamics is illustrated by showing the time evolution of the probability of 
\Q , occupation of the electronic levels for simple initial preparations of the state of the system. 
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Q ' PACS numbers: 42.50.-p, 32.80.Qk, 32.80.Wr 
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The Jaynes-Cummings model (JCM) Q is the most elementary quantum model for the interaction of matter with 
an electromagnetic field. In this model, a two-level atom is coupled to a single quantized mode of a cavity field 
and it is one of the few fully quantum-mechanical models that is exactly solvable. It exhibits some unexpected 
nonclassical features as the revivals of the Rabi oscillations in the atomic inversion, for instance 0. The quantum 
£SJ , origin of this revival is a direct consequence of the non-vanishing commutation relation between the creation and 
annihilation operator of photons in the field mode ,'2]. During the past decades, the JCM has been extended to more 
general Hamiltonians including multilevel atoms U multimode or external fields Q| , multi-atom configurations [f| , 
and multi-photon transitions @, just to mention a few examples. All these generalizations are part of what is called 
cavity quantum electrodynamics (cavity QED) and the quantized electromagnetic field plays a fundamental role in 
I . those settings. By the other hand, trapped ions interacting with classical fields have gained considerable interest 
qq ' in the past few years, mainly because the system dynamics closely resembles that of the JCM, with the quantized 
harmonic motion of the ion's center-of-mass playing the role of the field Q ■ Significant experimental advances in the 
^sO • generation of quantum states in such a system [g] have also indicated that trapped ions are very suited for the study 
and physical implementation of quantum dynamics Q . It seems clear that a system comprising quantum cavity fields 
■ and trapped ions undergoing quantized harmonic motion could bring many interesting consequences and potential 
fS |' applications. A typical setup would be that in which an ion trap is inserted in a high finesse optical resonator so that 
! 1 1 , the ion could interact with both quantized and classical external fields |1C( . Several authors have studied this new 
setting in the framework of single-photon transitions and also Raman transitions driven by the cavity field and 
a classical external laser However, we are not restricted to single-photon transitions in the realm of cavity 

| QED. Other processes, such as two-photon transitions, are particularly important and have potential applications, 
e.g., in the generation of non-classical states of the electromagnetic field, such as squeezed states ^3 or even photons 
with correlations that violate classical inequalities |l 5| . Nonlinear transitions in trapped ions interacting with classical 
laser fields have also been studied and applications suggested [iljj . 

This paper is concerned with the study of atom-field two-photon interactions when including harmonic atomic 
motion. We derive an effective Hamiltonian which not only describes electronic two-photon transitions but also 
contains a kind of phase-coupling between the vibrational center-of-mass motion, the electronic degree of freedom, 
and the cavity field. Another interesting feature of that Hamiltonian is the presence of Stark-shifts depending on 
the number operator of the harmonic motion. We then analyze the system dynamics by means of the probability of 
finding the ion in the electronic ground state, which is an easily accessible quantity in ion trap experiments We 
discuss how this quantity is influenced by the statistics of simple initial preparations of the quantum state of motion. 
In other words, we systematically compare our model with the important and extensively studied TPJCM without 
the motional effects 



I. THE MODEL AND RESULTS 



We consider a three-level ion in a cascade configuration interacting with a single-mode quantized electromagnetic 
field enclosed in a high finesse cavity. The schematic level structure is depicted in Fig.lJTJ. It is assumed two-photon 
resonance between the upper |e) and lower \g) atomic levels and the intermediate level |r) is kept off- resonant. The 
general Hamiltonian describing the interaction of trapped ions (or atoms) and electromagnetic fields is discussed in 
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several papers 0, El El 03 IH1 > an d in our case it reads 

i? = ^a^a + ojtfb + E e a ee + E r a rr + E g a gg + [gi(a gr P + a rg b) + g%{a re tf + a er b)] cos 77(0^ + a) (1) 

where a'(a) denotes the creation (annihilation) operator of the center-of-mass vibrational motion of the ion (frequency 
u), b'(b) is the creation (annihilation) operator of photons in the field mode (frequency cu), &ij = is a transition 
atomic operator, Ei is the energy of the atomic level |z), g\ e g 2 are the ion-field coupling constants for the transitions 
\g) — > I r) and \r) — ► |e), respectively, and r\ = 2i:a /X is the Lamb-Dickc (LD) parameter, being a the rms fluctuation 
of the ion's position in the lowest trap eigenstate, and A the wavelength of the cavity field. We have taken (h = 1) and 
this convention will be followed in the rest of our paper. The Lamb-Dickc parameter is the same for both transitions 
(coupling constants g\ and g 2 ) because both are driven by the same (cavity) field, i.e. with the same frequency (or 
k vector). The cascade level structure and the field frequency [see Fig.JTJ] have been chosen such that by making 
5 = E e — E r — u> = Eg — E r + u> 3> <7i , g% one can derive an effective two-photon Hamiltonian by means of the adiabatic 
elimination of the level |r). 




FIG. 1: Schematic diagram of the three-level ion interacting with a single mode quantized field with frequency lo c . 



The general Hamiltonian is highly nonlinear because the function cos r\(a) +a) contains powers of operators of the 
center-of-mass motion. Each term in the power series expansion will be dependent on powers of the LD parameter 
77. We can say generally that the higher the value of 77 the stronger will be the influence of nonlinear terms in the 
Hamiltonian QJ. It is well know that appropriate choices of the ion-field detuning 5 may lead to different kinds of 
couplings in the rotating wave approximation (RWA) 0- In general, it is assumed S = kv, with k integer, what leads 
to either transitions between the trap eigenstates (k-sideband Hamiltonian, for k ^ 0) or just a phase-coupling with 
no energy transitions for the ion's center-of-mass harmonic motion (carrier Hamiltonian, for k = 0). However, it is 
not necessary to have S = in order to forbid those transitions (make them unlikely). The same result applies for a 
less demanding situation in which 6 <C v. The advantage of having this less demanding condition on the detuning S 
is that now it is possible to think of a situation where the field and atom can be kept far off-resonance [6 3> 31,32) 
without having excitations of the center-of-mass motion. This may be achieved by setting v ^> S ^> gx, 92- In this 
case, it is possible to obtain a carrier two-photon Hamiltonian as we are going to show next. 
The Hamiltonian Q in the interaction picture reads 

Hi = [gi(a rg be~ ut + H.c.) + g 2 (a er be lSt + H.c.)} cos7 ? (a t e wt + ae~ 11 "), (2) 

which can be rewritten as 

Hi = Y}9i{°r g be- i[5+v[a - m + H.c.) + g 2 {a er be^ +vi - a -^ t + H.c.)]f(tf,% a, (3), (3) 

where 

f(tf,a;a,0) = ^r[(*vT +f3 + Hv) a+f V a ^. (4) 
2 aim 



3 



Analyzing the temporal dependence of Hamiltonian one can see that the frequencies may be carefully chosen 
allowing the RWA to be performed. In this approximation, only the slow frequency terms are kept while the rapidly 
oscillating ones are discarded. For sufficiently short interaction times, this approximation is quite accurate as far as 
the coupling constants are not too strong. How strong the coupling constants must be when compared to the other 
frequencies of the problem is usually found by using time dependent perturbation theory. In our problem, the terms 
in oscillate in time as e l± ( 5J rkv)t ^ w jth k integer. In the regime 6 -C v, the slowly oscillating terms are those with 
temporal dependence e ±ll5t . This only happens when a = (3 in @. Then, by dropping out the rapidly oscillating 
terms in (|3J), we arrive at the following approximate Hamiltonian in the original picture 

H = va'a + ujtfb + E e a ee + E r a rr + E g a gg + f{a^a)[{g\a gr + g 2 a re )b^ + (fliaVg + 92<r e r)b] (5) 

where, 

°° / i \Q„2a„ta„a 

/(fitft) = e-" 2 / 2 £ M) V / ° = e- 2 / 2 : J (2 V V¥&) :, (6) 
* — ' a'. 

and : Jo : is the normally ordered zeroth order Bessel function of the first kind. In an appendix in |l7j . it is shown 
how to convert : Jo : in a form that does not contain the normal ordering symbol but is an expansion in Fock basis. 
This is very useful for doing the plots presented in the last part of this paper, and the relation between both forms is 
in our case given by 

oo 

: J Q (2 V V¥&):=Y^L n (v 2 )\n)(n\, (7) 

71=0 

where L n is the Laguerre polynomial. The Hamiltonian (JSJ describes a situation in which the ion's center-of-mass 
motion couples to field and electronic degree of freedom in such a way that only phases are involved, i.e. there is no 
real transitions between vibrational energy levels. There is just photonic transitions taking place. The magnitude of 
this phase-coupling is contained in f(a'a) which is a function of the number operator of the vibrational motion. In 
the limit 77 — > the function /(oJd) tends to the identity operator, what corresponds to the free motion of the ion. 

In order to derive an effective two-photon Hamiltonian, let us now find the equations of motion for some relevant 
operators. The starting point would be a eg that is supposed to be present in the effective Hamiltonian once its 
function is to cause the ion to make a direct two-photon transition from the state \g) to \e). 

The Heisenberg equation for a eg using J5j) is given by 

* J^ e s = ( E a ~ E e)<r eg + 9if(a*a)a er tf - g 2 .f '(a 1 a)a rg b^ '. (8) 

The right hand side of JSJ) involves operators in the form f{c^a)aijW. We need to compute the Heisenberg equations 
for these operators as well, and they are given by 

i^[f(tfa)a er &] = (E r -E e -Lu c )f(a^a)a er P + [f(a^a)} 2 [ gi b r ba eg + g 2 b^(a ee - a rr )], (9) 

ij t [f{a)a)a rg P] = (E g - E r - uj c )f{a) a)a rg fr + [/(a 1 a)} 2 { gi tf 2 (a rr - a gg ) - g 2 (l + tfb)a eg ]. (10) 

The adiabatic elimination of the level |r) follows from (j^J) and l|10|) by considering the condition 6 3> gi,g 2 - To 
make this clear, it is convenient to define new operators in an appropriate reference frame as &ij — > a^e 1 ^ 3 ^®*'*, 
b(t) — * be~ Wct e a(t) — ► ae~ tvt . In this new frame, we may rewrite ©, © and (|ll)fl as 

i^a eg = 9l f(tfa)a er fre m + g 2 f(tfa)& rg Pe- iSt , (11) 



dt 

d 
Jt 
d_ 

dt 

We then integrate l(T^|l and l(T^|) under the assumption that 5 ^> gi,g 2 - By substituting the result of the integrations 
in 111|) and also considering that the level \r) is not initially populated, we end up with the following equation of 
motion in the original Heisenberg picture 



i^[f(tfa)* er fr) = [f(tfa)] 2 [gM& eg +g 2 tf 2 (* ee -a rr )}e- i5t , (12) 
-[f{tfa)& rg P] = [f(^a)] 2 [ 9l P\a rr -a gg )-g 2 (l + ^b)a eg ]e m . (13) 



«^eg = f(^d) 
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-f 2 {a)d){a ee - a gg p 2 (14) 
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which may be obtained from the effective Hamiltonian 

H = Hq + H Stark + Hj, (15) 

being the free part given by 

Ho = va'a + u c Wb + u} c (a e e — figg), (16) 

the Stark shifts 

2 2 

Hstark = ^/ 2 (ata)(l + Pb)a ee + 9 -jrf{tfa)tfbct gg , (17) 

and the two-photon interaction term between ion and field given by 

Hi = ^f 2 (a^a)(a eg b 2 + a g( P 2 ). (18) 

It describes two-photon transitions between the levels |e) e \g) with a coupling constant that depends on energy of 
the center-of-mass motion via f 2 (a*a). 

An interesting feature of the Hamiltonian l|15|) is the dependence of the Stark shifts upon the motion, fact that is 
mathematically expressed by the presence of f 2 (a^a) in l|17f) . This is the first remarkable feature of our model that 
is not present in the ordinary TPJCM. This means that the vibrational degree of freedom will have an important 
influence on the dynamics of the system not only due to the statistics of the quantum state of motion and the phase- 
coupling with the rest of the system but also via the Stark-shifts. As we are going to see below, even small variations 
of the Lamb-Dicke parameter will produce significant changes in the dynamics of the system. This makes the model 
here presented useful for the investigation of quantum aspects of the light-matter interaction. 

Once we have derived the phase-coupling two-photon Hamiltonian with motion-dependent Stark-shifts i|15|) , which 
is the main result of our paper, we now proceed with the study of the dynamics of the electronic levels for simple but 
interesting and experimentally feasible initial preparations of the system. In order to do that, we will compute the 
occupation of level \g) that is defined as 

P g (t) = \(g\m)\ 2 - (19) 

It follows from the form of l|15(l that as long as the ion is initially in its electronic ground state, the global state of 
the system may be written as 

oo oo 

\ip(t)) = ^2 y]am n (t)\m, n, e) +b mn (t)\m,n-2,g), (20) 

m— n— 

where the index m is referred to the motion and n to the field. If one finds a mn (t) an b mn (t), the desired probability 
P g (t) is just the summation in m and n of the squared absolute value of b mn (t). In the interaction picture, the 
coefficients obey the following system of differential equations 

i-^a m n(T) = Xl £w(t) + X2 bmn(T), 

i^b m n(T) = X3 b mn (r) + \2 a mn (r), (21) 
at 

where we defined 

_ .9152 _ , 9i _ P(m) 

9 = — r~ ; T = gt; r = — ; Xl = (" + 1)5 

o 92 r 

X 2 = V(n + l)(n + 2)f 2 (m); X 3 = rf 2 (m)(n + 2). 

(22) 

The solution of (|21J) gives one all the information available about the physical system described by the Hamiltonian 
l)15[l. For the purposes of this paper though, we are interested in the simplest case gx = g<i (r = 1). Also, we will be 
particularly interested in initial preparations given either by 

a C01 - e -l«l 2 /2_^ j 

Li mn yJ ) — c , Urip 

V TUl 

b mn (0) - (23) 
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t (0) = e -(l«l 2 +l0l 2 )/2 a P 



y/mln\ 

b mn {0) = 0. (24) 

In both cases the ion is initially in the electronic excited state |e), and particulary in (|23JI the motion is in the coherent 
state | a) and the field in the Fock state \p) containing exactly p photons, while in (|24fl the motion and field are initially 
prepared in coherent states \a) and \(3), respectively. The experimental generation of vacuum Fock and coherent states 
of motion for trapped ions have already been reported @ as well as of electromagnetic cavity fields 01 • The solution 
of l|21|) for r = 1 and the ion initially in the excited state is 



i(t) = a mn (0) 



f 2 (m) 

cos(A to „t/2) + i— sin(A m „r/2) 



b mn (r) = -2ia mn (0)- sin(A mn r/2), (25) 



where A mn = \// 4 (m) + 4%|. The sought probability is finally found to be 

oo oo a 2 

P M = E E K»(°)| 2 T^ sin 2 (A m „ r/2). (26) 

m=0 m=0 mn 

In what follows we will be studying the effect of the motion on P g (r) . Mathematically, the ordinary TP JCM may be 
obtained by making 77 — ► 0. We note that both the frequencies and amplitudes of each term in equation l|26(l have 
distinct contributions from the quantized field and from the ion's vibrational motion. For the first initial preparation 
(|23|l . we clearly see [Fig. 10] how strong the influence of the atomic motion is because slight changes in the Lamb-Dickc 
parameter 77 lead to quite different behaviors. The expected Rabi oscillations found in the TP JCM are modified as 
the parameter 77 is increased. For small values of 77, i.e. in the limit of the TPJCM, the coherent state of motion 
induces almost complete periodic collapses and revivals. By increasing 77, the dynamics changes again and looses 
that periodicity until a complete irregular pattern takes place. Similar strong effects of the harmonic motion on the 
field dynamics in cavity QED setups have been already reported, see for instance Di Fidio et al [l3[ . They consider 
a trapped ion in a Raman configuration interacting with a cavity field and an external laser. The main difference 
between our model and the one treated in [f3 | is that the two-photon configuration is in essence a periodic model 
while the Raman coupling does not necessarily lead to a periodic behaviour. Therefore, the inclusion of the atomic 
motion might not produce the same effects on those two different cavity-QED setups. In fact, the influence of the 
atomic motion should become more evident in the model treated here, as the oscillations induced by the atomic 
motion are in general not periodic, in contrast to the oscillations due to the interaction with the field. Considering 
now the second initial preparation l|24|) . for an initial coherent state for the field, the dynamics is also modified by the 
harmonic motion [Fig. 10]. In this case, the characteristic (almost) periodic evolution found in the in the ordinary 
TPJCM and reported in several papers Q, is again modified by the increasing of 77. The beats due to the statistics of 
the initial state of the center-of-mass motion (coherent state) clearly destroy the regular patterns. Moreover, for this 
initial preparation in which both the field and the center-of-mass motion of the ion are prepared in coherent states, 
we found for intermediate values of the LD parameter the interesting behavior of revivals ocurring at longer times 
|l9j , as shown in Fig. (0 . This super revival or revival of revivals is a revival at long times of the Rabi oscillations and 
the ordinary short-time revivals. That is another special feature happening for this model that is not present in the 
TPJCM. Once the dynamics in the TPJCM is almost periodic for short and long interaction times, the concept of 
super revivals has no meaning in this case. All these different effects indicate some of the several interesting processes 
that might arise when considering cavity quantum electrodynamics with trapped ions, specially the model proposed 
here. We note, from the examples above, that the effects of the quantized field and the atomic motion are somehow 
superimposed and have peculiar features, such as, for instance, a dependence of the Stark shifts terms on the atomic 
motion. 



II. CONCLUSIONS 



We have investigated the interaction of a trapped ion with the quantized cavity field via two photon transitions. 
Particularly in this paper, we have studied this system focusing on the fundamental aspects of the understanding 
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of light-matter interaction. We have derived an effective Hamiltonian under the rotating wave approximation and 
presented its analytical solution. This effective Hamiltonian contains motion-dependent Stark-shifts and a coupling 
constant which is a function of the intensity of motion, i.e. it is a function of the massive oscillator number operator. 
Both features are not present in the original two-photon Hamiltonian of the TPJCM. In particular, we have calculated 
the evolution of the population of the electronic levels, and found how the center-of-mass motion decisively changes the 
dynamics. In general, a treatment of the open system by including cavity losses, spontaneous electronic emission, and 
so forth, is necessary, although it rarely possesses a closed analytical solution. Some important and exceptional cases 
that admit analytical solutions for especial regimes were pointed out and treated by Di Fidio et al 01 ■ Even though 
the system considered in 0] is not identical to ours as they treat the case of Raman transitions, most their findings 
may be applied to our problem. As they point out, the inclusion of these realistic conditions will certainly have a 
destructive effect as long as quantum coherences are involved. In practise, some behaviors of P g (r) presented here will 
not be observed in the current experimental setups, and it is the case for the super revival phenomenon. However, it 
is always important to investigate the presence or not of these typical quantum effects. On the other hand, differences 
between our model and the ordinary TPJCM arise also at shorter times as, indicated in Fig[3and Figdl Therefore, 
it would be possible to observe some of those effects with the achievement of the strong coupling regime; although 
losses and dephasing effects are present, they will have a much weaker magnitude than the cavity-ion coupling itself. 
The system comprising cavities and trapped ions studied here is now under intense experimental investigation with 
significant control improvements pp| |. and efforts are being made in order to reach the strong coupling regime. 
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and FAPESP (Fundagao de Amparo a Pesquisa do Estado de Sao Paulo) grant number 02/02715-2, Brazil. 
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FIG. 2: Time evolution of electronic ground state population for different values of the Lamb-Dicke parameter. The system 
has initially been prepared in the state \tp(0)) = |e, a = 2,p = 0). 
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FIG. 4: Long time evolution of the electronic ground state population. The system has initially been prepared in the state 
|V(0)> = |e,a = 2,/3 = 2>. 



